Theoretical analysis of magnetic force microscopy contrast in multidomain states of 
magnetic superlattices with perpendicular anisotropy 
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Recently synthesized magnetic multilayers with strong perpendicular anisotropy exhibit unique 
magnetic properties including the formation of specific multidomain states. In particular, antiferro- 
magnetically coupled multilayers own rich phase diagrams that include various multidomain ground 
states. Analytical equations have been derived for the stray-field components of these multidomain 
states in perpendicular multilayer systems. In particular, closed expressions for stray fields in the 
case of ferromagnetic and antiferromagnetic stripes are presented. The theoretical approach provides 
a basis for the analysis of magnetic force microscopy (MFM) images from this novel class of nano- 
magnetic systems. Peculiarities of the MFM contrast have been calculated for realistic tip models. 
These characteristic features in the MFM signals can be employed for the investigations of the dif- 
ferent multidomain modes. The obtained results are applied for the analysis of multidomain modes 
that have been reported earlier in the literature from experiments on [Co/Cr]Ru superlattices. 
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PACS numbers: 75.70.Cn, 75.50.Ee, 75.30.Kz, 85.70.Li 



Magnetic multilayers with strong perpendicular 
anisotropy are currently investigated as crucial elements 
in magnetic sensors, storage technologies, and magnetic 
random access memory systems A large group of 
them belongs to systems with antiferromagnetic coupling 
through non-ferromagnetic interlayers (e. g. Co/Ru, 
Co/Ir, [Co/Pt]Ru, [Co/Pt]NiO superlattices) [10, 111- 
These nanoscale synthetic antifcrromagnets are charac- 
terized by new types of multidomain states, unusual de- 
magnetization processes and other specific phenomena 
@; S HI- In contrast to other bulk and nanomagnetic 
systems, the multidomain states in perpendicular an- 
tiferromagnetic multilayers are determined by a strong 
competition between the antiferromagnetic interlayer ex- 
change and magnetostatic couplings 0, H, 0- The re- 
markable role of stray-field effects in synthetic antiferro- 
magnets and the peculiarities of their multidomain states 
are currently investigated by high resolution magnetic 
force microscopy (MFM) (for recent examples of success- 
ful experimental tests on domain theory by MFM see, 
e.g. Refs. [1, [||]). From the theoretical side, only few 
results have been obtained on MFM images in antifcr- 
romagnetically coupled multilayers, mostly by numerical 
methods 0,^0- Here we present an analytical approach 
that provides a comprehensive description of stray-field 
distributions and MFM images in multidomain states of 
these nanostructures. We show that the stray-field com- 
ponents and their spatial derivatives, that are crucial for 
an analysis of MFM contrast, own distinctive features 
for different multidomain states. These features allow to 
recognize the particular distribution of the magnetiza- 
tion at the surfaces of domains and in the depth of the 
multilayers. The quantitative relations from theory for 
the MFM contrast can also serve to determine the values 
of magnetic interactions, i.e. materials parameters of an 
antiferromagnetic multilayer. We apply our results for 
an analysis of multidomain states observed in [Co/Pt]Ru 
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FIG. 1: (Color online) Equilibrium period of stripe states in 
multilayers. For various values of multilayer repeat numbers 
N, the stripe period D in units of the characteristic lengths I 
is shown for the ferromagnetic mode (dotted line), and anti- 
ferromagnetic mode (solid lines) as function of the magnetic 
layer thickness h/l for the ratio of interlayer to ferromagnetic 
layer thickness s/h = 0.1. The dashed line indicates the solu- 
tions for a single layer. The inset introduces the geometrical 
parameters of the model. 



multilayers @. 

As a model we consider strong stripes, i.e. so-called 
band domains in a superlattice composed of N identical 
layers of thickness h separated by spacers of thickness 
s, (see, Fig. [T]). Note, that the term "stripe domains" is 
also commonly used to denote multidomain patterns con- 
sisting of stripes with weakly undulating magnetization 
which, however, stays predominantly in the layer plane 
[Tol ]. On the other hand, the term band domains is used 
to describe structures of homogeneous domains with per- 
pendicular magnetization that alternates between up and 
down direction. These two types of stripe domains should 
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not be confused. In the model for an extended multilayer 
film, the multilayer is taken to be infinite in x— and y— 
directions. The stripes with alternating magnetization 
M along the z direction and with |M| = M = const 
have the period length D and are separated by domain 
walls of vanishingly small thickness (Fig. Q]). However, 
domain walls contribute a positive excess energy with an 
area density a, which is included in the model as one of 
the materials parameters. The ferromagnetic layers are 
antiferromagnetically coupled via a non-ferromagnetic 
spacer. This interaction imposes antiparallel orientation 
of the magnetic moments in adjacent layers, while mag- 
netostatic forces favor parallel orientation of the magne- 
tization. As a results of this competition three different 
ground states can be realized depending on the mate- 
rials and geometrical parameters of the multilayer Q. 
Namely, the homogeneous antiferromagnetic state, and 
stripe domains with parallel or antiparallel arrangement 
of the magnetization in adjacent layers (Fig. [5]). We de- 
note these latter two multidomain modes as ferro and 
antiferro stripes. The energy density for both types of 
stripe phases can be written as 



The first term in wn is the domain wall energy, J > is 
the antiferromagnetic exchange interaction via the spacer 
layer. The stray field e nerg y w m can be expressed as a 
set of finite integrals 0, El], 

1 rT ~ 1 

w m (D) = 1 + fi(£>, h) T -T7 E ( N - fc ) E »( D > h > s )' ( 2 ) 



k=l 



where T = h + s is the period length of the superlattice, 
S fe = 2 fl(D, Tk) - fl(D, Tk + h) - fl(D, Tk - h), and 
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The upper (lower) sign in Eqs. (H|), and in @ for w m , 
corresponds to ferro(antiferro) stripes. The equilibrium 
domain configuration of the stripes is derived by mini- 
mization of wn with respect to the stripe period D Q. 
Introducing a new length scale, based on the character- 
istic length I = a/ (AttM 2 ), one can express the solutions 
for reduced periods D/l as functions of three parameters, 
the reduced geometrical sizes (h/l), (s/l), and the repeat 
number of the multilayer N. Typical solutions D/l as 
functions of h/l are presented in Fig. [TJ The solutions 
for ferro stripes exist for arbitrary values of h and s. In 
the two limiting cases of large and small spacer thickness, 
the solutions asymptotically approach the behavior of the 
known solutions for individual layers (see Ref. [12, [Hj]) 
with a thickness h for the case s 3> h and with effec- 
tive thickness hN for s -C h. The solutions for antiferro 
stripes with even N exist only in an interval bounded 
from below, h > h cr (N). The period D tends to infinity 
at a critical thickness h cr . Strictly speaking, for odd N 
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FIG. 2: (Color online) Calculated stray-field profiles Hz (x) 
for ferromagnetic (a) and antiferromagnetic (b) modes for 
[[Co/Pt]7CoRu]4 multilayers investigated in c|. Insets show 
the perpendicular stray-field component H^ n \zo) at the cen- 
ter of the stripes in dependency on the distance zo above the 
multilayer surface. In ferro stripes this function monotoni- 
cally decreases with increasing z, while in antiferro stripes it 
has a maximum at a finite distance from the surface. 



antiferro stripes (similarly to ferro stripes) exist for ar- 
bitrarily small layer thickness h. However, their periods 
increases so steeply that a single domain state is prac- 
tically reached when the period exceeds the lateral size 
of the layer. The calculated periods for odd-numbered 
multilayers has similar have similar size as those for an- 
tiferro stripes with even N. The phase diagrams con- 
sidering these stripe ground-states show that both types 
of stripes can exist as stable or metastable states in ex- 
tended and overlapping ranges of the material parame- 
ters 0j- The extended co-existence regions of different 
types of multidomain states in the phase diagrams also 
entails the possibility to create complex "interspersed" 
patterns that consist of subdomains with ferro and anti- 
ferro stripes. For identical values of the materials param- 
eters the equilibrium domain widths for ferro and anti- 
ferro stripes can differ considerably (see Fig. [T|). Hence, 
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the mixed stripe patterns can include regions with dif- 
ferent domain sizes. Exactly such structures have been 
observed in some [Co/Pt]Ru multilayers [14]. In addi- 
tion to the differing characteristic periods of ferro and 
antifcrro stripes, these stripe patterns also cause differ- 
ent distributions of the stray fields H 1 ™) at the sample 
surfaces. The stray-fields can be probed by magnetic 
force microscopy, however, the properties of the stray 
fields peculiar to the different stripe patterns are rather 
subtle. Thus, the experimental observation and quanti- 
tative evaluation of these differences must be based on a 
detailed comparison with theoretical model calculations. 

Convenient analytical expressions for stray field com- 
ponents and their spatial derivatives for multilayers with 
ferro and antiferro stripes can be derived from Eqs. {T])- 
@. The mathematical equivalence between electro- and 
magnetostatic problems allows to treat a magnetized 
layer as a system of "charges" distributed over its surfaces 
[lfj ] . The stray field above the sample surface can be ex- 
pressed as a superposition of the stray fields from the 2N 
interface planes with "charged" stripes (Inset in Fig. [1]). 
By solving the magnetostatic problem for a plane with 
"charged" stripes (see Appendix) one can write the fol- 
lowing solutions for the stray field components h^™ 1 ' [x, z) 



where 



^ m) (a;,z) = 2Mln 



cosh (2nz/D) - sin (2ttx/D) 



cosh (2ttz/D) + sin (2nx/D) 



, (4) 



h{ m) (x,z) = 4Marctan [cos (2nx/D) /sinh (2-kz/D)] .(5) 

v (x,z) 

Then the total stray field of the multilayer at a distance 
Zq above the surface can be written as 

N—l 

U^(x,z ) = r * \h {m) (x,z +Tk) 



-h (m) (x,z + Tfc + /i) 



(6) 



The factor Tfc = (— l) fe holds for antiferro stripes, and 
Tk = 1 for ferro stripes. Spatial derivatives of Hi' with 
respect to z are important for the analysis of the MFM 
images. The derivative T„(x, zq z™can be 

derived analytically by differentiation of Eq. © as 

2vr 



T n (x,z )=4M^ — 

N-l 

x X! ^ k \. v ™( x > z ° + Tk ) ~ v n(x, z +Tk + h)] . (7) 

We introduce here a set of functions v n (x, z) which are 
derivatives of the function vo(x,z) defined in Eq. (|5|) 
with respect to the normalized geometry parameter £ = 
2-kz/D, 

v n (x, z) = d n v /dC = cos {2-kx/D) (8) 
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[cosh {Attz/D) ± cos (4irx/D)} /2, 

- cosh (2irz/D) , 
sinh (2ttz/D) (1 + .g_), 

- cosh(2irz/D) (2g 2 _ - g\ + 2g+ - 2), 
sinh (2nz/D) 

5.9-)] 



x [6(. 9 _ + l) 2 (. 9 _ -1) + <?+(! 



(9) 



Together with the equation dwiy/dD — 0, which de- 
termines the equilibrium domain period, Eqs. ([6]) and 
([7]) describe the stray field and its spatial derivatives as a 
function of the coordinates x, z for a multilayer in a stripe 
state. The stray field H^ m )(a;, zq) ^ and the derivatives 
T n (x, zo) ([7]) are expressed as sets of analytical functions 
Q, © ©, and ©. The expressions depend on the 
geometrical parameters and, via the equilibrium domain 
widths D, on the material parameters of the multilayer. 
These analytical expressions can be readily evaluated by 
elementary mathematical means. 

In order to demonstrate the main features of the 
stray fields Hi m \x,z) from the ferro and antiferro 
stripes we evaluate these functions for a multilayer 
[[Co/Pt]7CoRu)]4 with magnetic and geometrical param- 
eters corresponding to a sample that was investigated 
experimentally in [3(. In this superlattice the ferromag- 
netic constituents are magnetic [Co/Pt]7 multilayers with 
thicknesses of the ferromagnetic Co-layers 0.4 nm and 
thickness of the Pt layer 0.7 nm. The non-ferromagnetic 
Ru spacer has a thickness of 0.9 nm and mediates an indi- 
rect antiferromagnetic interlayer exchange. The domain 
period has been determined as D — 260 nm [3] . For ferro 

and antiferro stripe modes the functions Hz {x, z) are 
markedly different both in the intensity and in the loca- 
tion of characteristic extremal points (Fig. [2] ) . Moreover 
they display qualitatively different functional dependen- 
cies on the distance from the multilayer surface zq (see 
Insets in Figf5]). 

The stray-field distribution over the multilayer surface 
can be viewed as a superposition of magnetostatic fields 
from systems of "charged" bands. This allows to give 
a simple physical interpretation for the main features of 
the stray field profiles in Fig. [21 First we consider ferro 
stripes. It is convenient to separate the total stray-field 
over a domain to two contributions: one created by the 
top and bottom bands of the domain(" self ' field), and 
that produced by all other "charged" bands. In the do- 
main centers near the sample surface, zq <C D, the stray 
"self" field is small due to screening effects of the domain 
top and bottom surfaces. Because the bands change their 
"charges" at the domain walls the stray field is substan- 
tially enhanced above the walls. As a result, the profiles 
ffl" 1 ' (x, z) have characteristic wells in the domain centers 
for zq ^ D. For increasing distance zq from the surface 
the difference between values of Hi" 1 ^ in the center and 
at the domain edges decreases due to the increasing in- 
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fluence of neighbouring poles. For large distances zq the 
wells disappear and the profiles obtain a typical bell-like 
shape (compare the traces for zq = 10 and 30 nm in 
Fig. [2] (a)). The antiferro stripe mode can be obtained 
from those for ferro stripes by changing the magneto- 
static "charges" for the bands in all even layers. This 
weakens the total stray field and sharpens the difference 
between the stray fields at the center and near the do- 
main edges. Finally, the competing character of the stray 
field contributions from odd and even layers causes the 
nonmonotonic dependence of H z m \z) (Inset in Fig. [2] 
(b))- 

Generally the functions (x , zq) , 



t! F) (:e, Zq), 

(x, zq) have a number of characteristic features that 
can be utilized in new methods to investigate the mul- 
tidomain modes. One method can be based on mea- 
suring the MFM contrast in the center of the domains. 
In this case x — kD, k — 0,1,2,3... and the functions 
v(z) = v n {kD,z) are reduced to the following expres- 
sions 

Vo(z) — arccos [tanh (2nz/D)] , 

vi(z) = - cosh -1 {2-kz/D) , 

v 2 {z) = sinh {2-kz/D) cosh -2 (2-kz/D) , 

-3 



v 3 {z) = -(1 - sinli {2-kz/D)) cosh" J (2ttz/D) 



V4,(z) = sinh (2irz/D) 



sinh 2 (2nz/D) - 
cosh 4 (2ttz/D) 



(10) 



Profiles v n (zo) from Eqs. {TDJ for n = 1,2,3,4 are plotted 
in Fig. [3] Their characteristic features can be used to as- 
certain the type of the stripe mode and should even allow 
quantitative evaluation of magnetic properties in multi- 
layers, if and when they are accessible in experiment. 

Owing to the properties of the magnetic probe, a sig- 
nal from a magnetic force microscope generally differs 
strongly from the H z (x) profile. In order to compare the 
expected outcome of different domain configurations, the 
MFM contrast has to be calculated for realistic tip mod- 
els. The MFM signal for a magnetic cantilever oscillating 
in z-direction is given by [la, [16| 



dz 



T M0 



/ ip 



d 2 [M(* l P)(r) -H( m )(r) 
dz 2 



dV.{ll) 



Here, A$ is the measured phase shift between excitation 
and oscillation due to the force gradient dF z /dz that 
acts on the cantilever in the stray field of the sample 
Hr")(r). Q and k are the quality factor of the oscil- 
lation and the spring constant, respectively. Assuming 
a rigid tip magnetized in z-direction, i.e., a tip with a 
homogeneous magnetization distribution, M^ tip ' = const 
throughout the tip volume, that does not change during 
the scan across the stray field of the sample, the expres- 
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FIG. 3: (Color online) Characteristic functions v n {zo) = 
Vn(0, Zq) (Eq. (flu) ') describe MFM images in the center of 
domains. Inset indicates the location of the extrema of the 
functions v n (x, zq) in the xOz plane. 



sion for the force gradient simplifies to 



dz 2 



/ ip 



d 2 H z m \r) 
dz 2 



dV. 



(12) 



The volume integration is a crucial step as it modifies the 
signal compared to the profile estimated by the second 
stray field derivative. 

A realistic tip geometry can be modelled by a trun- 
cated triangle placed in the a:-z-plane (see inset in Fig. [6]). 
This mimics the two parallel sides of the 4-sided pyrami- 
dal geometry of a typical MFM tip. The two-dimensional 
model simplifies considerably the calculations. The error 
incurred by this reduced model is minor because of the in- 
finite extension of the domain models in the y-direction. 

To demonstrate this approach we analyze the MFM 
contrast measured across a thick [(Co/Pt)g)CoRu]i8 mul- 
tilayer prepared at Hitachi GST (for details on these mul- 
tilayers, see Refs. [3j). The MFM pictures show a typical 
maze pattern of ferro stripes with perpendicular magne- 
tization at room temperature Fig. [?] (a)). For compari- 
son with contrast calculations the MFM signal along the 
marked line of the image was recorded repeatedly (with 
the slow scan axis disabled) at lift heights of 10, 20, 30, 
40, 50 and 60 nm. Averaged scan lines for each lift height 
in Fig.[4](b) display a reduced contrast for increasing scan 
height. 

The very regular domain pattern in the center of this 
scan (framed area) was modeled according to Eq. (J7J) as 
parallel FM stripes with the period of 360 nm based on 
an N = 18 multilayer with the known layer architec- 
ture. With the resulting second stray field derivatives the 
MFM phase shift was computed according to Eq. us- 
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x (|im) 



FIG. 4: (Color online) Maze domain pattern in 

[(Co/Pt)8)CoRu]i8 observed by MFM at room temperature 
(a) and line scans along the marked line for zo = 10, 20, 30, 
40 ,50 and 60 nm. 



ing the above mentioned tip model. As cantilever param- 
eters a spring constant of k — 2 N/m and a quality factor 
of Q = 100 were used. The height, pyramid angle and tip 
apex were chosen as 4 /*m, 30 and 50 nm, respectively, 
and the film coating was assumed to be 30 nm. As an ad- 
justable parameter the tip magnetization Mz, tip was set 
to 3T0 6 A/m. With this reasonable value the computed 
phase shift magnitude, line profile and lift height depen- 
dency compare very well with the experimental data (Fig. 
[5]). Such calculations can thus be used to predict differ- 
ences in the MFM contrast of the two distinct types of 
stripe domains. 

Fig. [5] shows the calculated force gradient profiles for 
the ferro and antiferro stripes presented in Fig. O The 
profiles reflect general features of the stray-field distribu- 
tions in the multidomain patterns. Here as well, quan- 
titative and qualitative differences allow to distinguish 
ferro and antiferro stripes. The signal from the anti- 
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FIG. 5: (Color online) Experimental (points) and calcu- 
lated (solid lines) force gradient profiles [(Co/Pt)8)CoRu]ig 
for MFM scans corresponding to the framed area in Fig. |4f b) . 



ferro stripes is clearly weakened, but it is large enough 
to be measured in a typical MFM setup which allows the 
detection of a few 10 /xN/m (10~ 2 dyn/cm) (T3| . How- 
ever, as quantitative MFM measurements are still rare 
and require precise calibration routines [16|, [17| the abso- 
lute value of the signal is not a reliable criterion for the 
distinction between different stripe configurations. More 
importantly, the force gradient directly above the cen- 
ter of a domain shows a monotonically decreasing sig- 
nal strength for increasing scan height zq in case of the 
ferro stripes. Above the antiferro stripes, on the other 
hand, the force gradient is increasing with increasing scan 
height, at least in the range from 10 to 30 nm. This qual- 
itative difference is demonstrated again in Fig. [5J where 
the force gradient experienced by a realistic tip model 
is plotted as a function of zq for the two cases. The 
non-monotonic behavior of the force gradients observed 
above a multidomain structure can be taken as a clear 
fingerprint of an antiferro stripe state within a multi- 
layer stack. The calculation even reveals a sign change 
as an additional signature of the antiferro stripes, but 
this appears at a scan height smaller than 10 nm, which 
is experimentally very difficult to access, as topographic 
information can influence the measurement. 

The solutions for the stray field ^ and its derivatives 
([7]) can be simplified in the practically important case 
of large domains, D ^> h. Expansion of T n (x,zo) with 
respect to the small parameter h/D -C 1 yields for ferro 
stripes 

T$P(x,z) = -4MN(2ir/D) n v n+1 (x,z), (13) 

and for antiferro stripes with even N 

T^ F \x,z) = -2MN{2Tr/D) n (l + s/h)v n+2 (x,z).(U) 

For antiferro stripes in multilayers with odd N the func- 
tions T n (x, z) are given by Eq. (fT3")) . Note that the func- 
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FIG. 6: (Color online) Calculated force gradient profiles for 
a 4-sided pyramidal tip and ferro (a), antiferro (b) modes 
in the center of domains as function of the distance zo from 
the surface (c) for [[Co/Pt]7CoRu]4 multilayers with the same 
geometrical parameters as in Fig. 2. 



tions Tn from the Eq. ([13)) are proportional to v n+ i, 
while the functions T^ 4 ^ in Eq. (jT5J) can be expressed 
by derivative of the functions u n+1 through the relation 
i>n+2 = dv n+ i/ dz. This means that in this limit of large 
domains, the functions T„ in the antiferro mode behave 
as z-derivatives of the corresponding functions of ferro 
modes. In particular for n = 0, the Eqs. (| 13[) and (fl~4"l) 
give the perpendicular stray field components for ferro 
stripes and antiferro stripes with odd N, correspondingly, 
by the expressions 

Hf\x,z ) = —AMN v x {x,z a l 

Hi AF \x, z ) = -2MN{1 + s/h) v 2 {x, z ). (15) 

In conclusion, we have presented analytical solutions 
for the stray field (Eq. ([6])) and its spatial derivatives 
(Eq. frj)) in multidomain states of magnetic multilayers 
with out-of-plane magnetization. These solutions can be 
applied to calculate the MFM contrast for a realistic tip 
geometry using Eq. (TP2|) . It is shown that the ground- 
state ferro and antiferro stripe structures in antiferro- 
magnetically coupled multilayers differ by their period 
lengths and by the spatial distribution of their stray- 
fields. Our analytical calculations executed within a sim- 
plified model of one-dimensional multidomain patterns 
with fixed magnetization orientation and infinitely thin 
domain walls are able to reproduce the general features 
of MFM images from antiferromagnetically coupled mul- 
tilayers. These features can be used to identify different 
types of multidomain patterns and extract values of the 
magnetic interactions. They also create a basis for more 
detailed investigations on more realistic models. Such 
models should consider distortions of the magnetization 
with tilting away from the perpendicular direction in size- 
able fractions of the domains due to the finite strength 
of uniaxial anisotropy, or a finite width of the domain 
walls between domains combined with the appearance of 
magnetic charge distributions at the walls. Additionally, 
it was recently shown that ferro stripes are unstable with 
respect to a lateral shift of domains in adjacent layers [l8j. 
This creates specific multidomain patterns with distinct 
modulations across the multilayer stack. Such effects 
can be consistently taken into account by numerical solu- 
tions of the micromagnetics equations for specific cases. 
Here, we have analyzed the backbone structure for all of 
these types of spatially inhomogeneous distributions of 
the magnetization within the multilayer and stray fields 
over its surface. 
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Appendix 

The scalar potential of a sheet with "charged" stripes 
(Inset Fig. [T]) can be derived by solving Poisson's equa- 
tion (see, e.g., [II, [l3|]) 



(x,z) 



AMD 



Esin(7m/2) . , , , , , 

— v 2 ' ' cos(2imx/D)exp(-2nnz/D) ■ (16) 



Using the identity n~ m = [(m - l)!]" 1 ^^e'^dt 



the Eq. (fl6l) can be transformed into the following form 



8MD f cos(2irx/D) cosh(t + 2ir z/D)tdt 



(x,z) 



7T J coah. z (t + 2nz/D) -sin 2 (2irx/D) 



(17) 



From this closed expression, the components of the stray 
field h.( m ) = — V</> are readily derived in the analytical 
form given by Eqs. (|4]) and ([5|). 
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